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摘 要 : 为 进一步 探讨 降水 量 和 水 土 保持 措施 对 径流 泥 沙 的 影响 ,深入 了 解 黄河 中 小 河流 水 沙 变 化 驱动 因素 , 采 


数 和 输 沙 模 数 呈 显 著 减 少 趋势 (P<0.05) ,水 土 保持 措施 


结果 表 昌 


了 M-K 检 验 .Morlet 小 波 线性 回归 和 结构 方程 模型 等 方法 ,对 1957 一 2021 年 关 川 河上 游 降 水 量 和 水 土 保 持 措施 面 
积 变化 及 其 对 径流 泥 沙 的 作用 机 理 和 过 程 进 行 了 研究 。 


:年 降水 量 呈 不 显著 减少 趋势 (P>0.05) ,径流 模 


下 积 呈 增 加 趋势 。 径 流 模 数 . 输 沙 模 数 对 降水 量 的 响应 由 


强 转 弱 , 对 水 土 保持 措施 的 响应 逐渐 增强 ;降水 量 和 水 土 保持 措施 影响 水 沙 变化 的 关键 路 径 是 "水 土 保持 措施 面积 
一 径流 模 数 ”, 其 次 是 “水 土 保持 措施 面积 一 径流 一 输 沙 模 数 ”; 再 次 是 “降水 量 一 径流 模 数 ”, 最 后 是 “降水 量 一 径流 
模 数 一 输 沙 模 数 ”; 降 水 量 对 水 沙 变 化 的 总 效应 分 别 为 0.42 和 0.38, 水 土 保持 措施 面积 对 水 沙 变 化 的 总 效应 则 分 别 
为 -0.72 和 -0.65 ,影响 流域 水 沙 变 化 的 主要 因素 是 水 土 保持 措施 面积 ,次 要 因素 为 降水 量 。 

关键 词 : 水 沙 变化 ; 水 土 保持 措施 ;结构 方程 模型 ， 关 川 河 上 游 


径流 量 , 输 沙 量 是 河川 径流 重要 的 水 文 要 素 ， 
是 评估 区 域 水 文 环境 动态 变化 的 重要 指标 ,也 是 流 
域 治 理 决 策 的 主要 依据 。 随 着 全 球 范围 内 气候 变 
化 和 人 类 活动 的 加 剧 ,流域 水 文 过 程 受到 了 不 同 程 
度 的 影响 ,进而 影响 流域 输 沙 量 ”;。 黄 河 是 世界 上 
含 沙 量 最 大 的 河流 ,黄河 水 沙 关系 变化 影响 到 流域 
内 社会 ,经济 .生态 环境 等 众多 领域 ,是 学 术 界 关注 
AYP), Zheng 等 "通过 对 比 黄土 高 原 122 个 水 
文 站 径流 泥 沙 数据 发 现 2008 一 2016 年 平均 径流 量 
和 输 沙 量 相 比 1971 一 1987 年 分 别 减少 22% 和 74%， 
并 运用 归 因 分 析 方 法 认为 人 类 活动 的 贡献 率 >72%; 
Zhang 等 "运用 回归 分 析 方 法 研究 祖 历 河流 域 
1955—2013 年 的 水 沙 演 变 规律 ,结果 发 现 降 水 量 、 
径流 量 和 输 沙 量 呈 逐年 减少 趋势 ;Rustomji 等 … 通 
过 分 析 1950 一 2000 年 黄土 高 原 流域 的 水 沙 关系 ,发 
现在 高 径流 条 件 下 含 沙 量 受 水 土 保持 措施 影响 呈 
大 幅度 降低 ;Zhang 等 所 研究 发 现 ,不 同时 间 尺 度 的 
水 沙 关 系 受 生态 建设 等 水 土 保持 措施 的 直接 影响 ， 
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不 同 径流 的 合 沙 量 均 呈 显著 降低 趋势 ; 雇 建 华 等 
从 高 含 沙 水 流 等 角度 分 析 水 沙 行为 特征 ,结果 表明 
水 土 保持 措施 及 植被 改善 降低 了 高 含 沙 水 流 发 生 
频率 ; 赵 娟 "采用 M-K 和 PP 突变 点 检验 、 双 累积 曲线 
法 分 析 了 黄河 支流 佳 芦 河 、 秃 尾 河流 域 水 沙 关系 变 
化 特征 及 与 人 类 活动 的 影响 关系 ,结果 表明 人 类 活 
动 对 径流 量 和 输 沙 量 的 贡献 率 大 于 降水 量 对 其 的 
贡献 率 ; 称 兴 民 等 5 在 变量 共 线 性 分 析 基 础 上 提出 
了 流域 降水 量 标 度 和 水 土 保 持 标 度 及 其 计算 方法 ， 
并 建立 流域 降水 -水 土 保持 -径流 统计 模型 ,在 佳 芦 
河和 般 尾 河流 域 水 土 保持 使 流域 径流 量 平均 减少 
10%~22% , 佳 芦 河 降水 变化 及 水 土 保持 措施 影响 
分 别 为 25% 和 75% ,而 秃 尾 河 分 别 为 35% 和 65%， 
同时 穆 兴 民 等 "研究 表明 ,人 类 活动 是 20 世 纪 70 
年 代 特 别 是 进入 21 世纪 以 来 黄河 泥 沙 减少 的 主要 
原因 ,20 世 纪 黄 河 泥 沙 减 少 主 要 是 由 于 水 利水 土 保 
持 工程 的 拦蓄 作用 ,21 世纪 以 来 的 减少 主要 为 规模 
化 恢复 植被 的 坡 面 阻 控 作用 。 径 流 泥 沙 的 变化 受 
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多 种 因素 影响 且 过 程 复 杂 , 以 往 人 研究 多 集中 于 不 同 
因素 对 径流 泥 沙 的 直接 影响 分 析 , 对 其 作用 机 理 和 
路 径 分 析 则 较 少 。 结 构 方 程 模型 作为 一 种 验证 性 
多 元 统计 分 析 技 术 , 有 助 于 观测 降水 量 水 土 保持 
措施 面积 与 径流 输 沙 之 间 的 直接 影响 、 间 接 影响 、 
总 效应 及 路 径 分 析 。 为 此 ,本 文 在 以 往 的 研究 基础 
上 ,选取 关 川 河上 游 1957 一 2021 年 水 文 站 点 降水 
Ht .水土 保持 措施 面积 、 径 流量 和 输 沙 量 数据 ,从 流 
域 水 沙 变 化 特征 、 水 沙 变化 归 因 分 析 及 水 沙 变化 对 
降水 和 水 土 保持 措施 的 响应 过 程 (路 径 ) 三 个 方面 
进行 研究 和 探讨 ”2 。 


1 数据 来 源 与 研究 方法 


1.1 研究 区 概况 

关 川 河 属于 祖 厉 河 的 一 级 支流 ,位 于 甘肃 省 中 
部 地 区 , 关 川 河上 游 包含 东 河 西河 两 条 支流 ,流域 
总 面积 1640.04 km’, 其 中 包含 东 河 流域 面积 791 
km ,西河 流域 面积 637 km ,两 河 在 安定 区 城关 镇 汇 
合流 入 关 川 河 主 河道 ; 东 河 与 西河 自 安 定 区 南边 边 
界 向 北 地 势 走低 ,侵蚀 地 表 形 成 黄土 丘陵 .堆积 河 
SAREE SF HS HY, JRE BR Ee YAK ; 气 
候 属 温带 季风 型 气候 ,气候 干燥 ,干旱 少雨 ,平均 降 
水 量 仅 413 mm 且 时 空 分 布 不 均 ,主要 集中 于 7 一 9 
H ,年 蒸发 量 1400 mm 以 上 ;土壤 主要 为 黄 绵 土 \. 灰 
钙 土 。 人 研究 区 以 水 力 侵蚀 为 主 ,1957 一 1980 年 年 平 
均 土 壤 侵 蚀 模 数 为 4558.8 t-km?-a!, BE 2021 年 
底 , 累 计 治 理 水 土 流失 总 面积 13.05x10: hm ,其 中 ， 
兴修 梯田 为 60963.01 hm ,水 保 林 为 37750 hm ,经济 
林 为 5380 hm ,人 工种 草 为 19835 hm? ,封禁 为 6560 
hm。 研 究 区 先后 设 有 网 口水 文 站 (1956 年 建站 一 
2000 年 )、 东 河 (定西 ) 水 文 站 (1984 年 建站 ) 和 西河 
水 文 站 (2001 年 建站 )3 个 水 文 观测 站 .6 个 雨量 观测 
站 。 水 文 站 雨量 站 点 分 布 如 图 1 所 示 。 
1.2 数据 来 源 及 预 处 理 

1957—2021 年 年 降水 量 径流 量 、. 输 沙 量 数据 
由 甘肃 省 水 文 水 资源 局 提供 。(1) 年 降水 量 , 人 研究 区 
共 设 有 东 河 西河、 内 官 营 红土 . 青 风 、 风 口 6 个 雨 
量 观测 站 ,平均 年 降水 量 采用 泰 森 多 边 形 法 计算 。 
(2) 径流 量 和 输 沙 量 数据 ,1957 一 2000 年 径流 量 和 
输 沙 量 数据 来 源 于 网 口水 文 站 观测 数据 ;2000 一 
2021 年 的 上 述 数据 来 源 于 东 河 西河 水 文 站 观测 数 
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图 1 关 川 河上 游 地 理 位 置 


Fig. 1 Geographical location of the upper Guanchuan River 


据 。 为 最 大 程度 消除 风口 水 文 站 迁移 带 来 的 影响 ， 
将 原 径 流量 与 输 沙 量 换算 为 年 径流 模 数 和 输 沙 模 
数 ,作为 基础 分 析 数 据 。 

水 土 保 持 措施 实 有 量 数据 来 源 于 甘肃 省 水 土 
保持 年 报 (1957 一 2021 年 ) 安定 区 水 土 保持 年 报 
(1957 一 2021 年 ) “甘肃 省 水 土 保持 综合 治理 措施 
效益 研究 (1979 一 2012 年 ) 29 等 科研 成 果 数 据 EC 
定 区 项 目 土地 详 查 二 调 .三 调 数据 .安定 区 林业 局 
详 查 数据 (2009 年 .2014 年 .2016 年 .2021 年 ) 及 退 
耕 还 林 工 程 (2000 一 2006 年 .2014 一 2021 年 ) 数 据 ， 
利用 上 述 资 料 修正 得 到 各 时 间断 面 项 目 区 实际 水 
土 保 持 措 施 面 积 。 根 据 不 同时 段 水 土 保持 措施 保 
存 率 ,计算 得 到 本 研究 区 历年 水 土 保持 措施 实际 
1.3 研究 方法 

(1) 采用 非 参 数 统计 检验 法 Mann-Kendall(M- 
K)22 分 析 关 川 河 降水 量 径流 模 数 与 得 沙 模 数 的 序 
列 变化 趋势 ,并 检验 其 变化 的 显著 性 。 

(2) 采用 Morlet 小 波 进行 小 波 分 析 ,根据 二 维 
等 值 线 闭合 中 心 的 正 负 值 和 小 波 系数 的 零点 来 判 
斯 降水 量 的 增 减 性 和 突变 性 ,其 主要 变化 周期 参照 
峰值 来 确定 。 

(3) 采用 线性 回归 分 析 降 水 量 .水土 保 持 措 施 
面积 与 径流 模 数 . 输 沙 模 数 的 关系 。 

(4) 采用 结构 方程 模型 (Structural Equation 
Modeling, SEM) 对 降水 量 水土 保 持 措 施 面 积 与 径 
流 模 数 输 沙 模 数 之 间 的 关系 进行 测度 和 分 析 , 建 
立 降水 量 和 水 土 保持 措施 对 径流 模 数 输 沙 模 数 关 
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系 式 ,最 后 确定 模型 参数 值 ””。 结 构 方程 模型 包 
含 结构 方程 和 测量 方程 。 结 构 方 程 用 来 测度 潜 变 
量 之 间 的 关系 ,测量 方程 用 来 测度 观测 变量 与 潜 变 
量 的 关系 。 公 式 (1) 和 公式 (2) 为 测量 模型 ,公式 
(3) 为 结构 模型 。 


X=Axf+6 (1) 
Y=Aynte (2) 
n=Bn+Té+t (3) 


式 中 :为 外 因 测 量变 量 向 量 ;Ax 为 在 it 上 的 因子 
WREE; 6 为 测量 误差 ;7 为 内 因 测量 变量 向 量 ;7 
为 内 因 潜 变 量 向量 ;4y 为 了 在 7 上 的 因子 荷载 矩阵 ; 
e 为 测量 误差 ;nm 为 内 因 潜 变量 向 量 ;é 为 外 因 潜 变量 
向 量 ;B 为 内 生 潜 变量 间 的 回归 系数 ;IT 为 外 因 潜 变 
量 对 内 因 湾 变量 之 间 的 回归 系数 ;为 测量 误差 。 为 
确定 模型 拟 合 指标 检验 模型 的 拟 合 优 度 ,选取 比较 
拟 合 标 准 化 估计 值 (Estimate) RUER (S. E.) . 近 
似 误差 方 根 (C. R./T) 和 P 值 4 个 指标 来 评价 模型 拟 
合 效果 所 。 模 型 的 构建 及 拟 合 检验 在 AMOS 21 和 
SPSS 26 中 进行 。 


2 结果 与 分 析 


2.1 水 沙 变 化 特征 

2.1.1 径流 模 数 年 际 变化 特征 “1957 一 2021 年 关 川 
河上 游 年 均 径 流 模 数 为 11466.71 mm .km a, BRK 
{HN 45221.97 m .km + a! (1967 年 ) ,最 小 值 为 
564.33 m°- km” - a`! (2011 年) ,最 大 年 径流 模 数 是 最 
小 年 径流 模 数 的 80.1 倍 。 由 径流 模 数 线性 拟 合 可 
知 ( 图 2a) ,年 径流 模 数 总 体 呈 显著 下 降 趋势 (P< 
0.05). 1957—1967 年 径流 模 数 呈 增 加 趋势 ;1967 一 
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2021 年 径流 模 数 逐年 减少 ,其 中 ,1999 一 2017 4F 48 
流 模 数 显著 下 降 (P<0.05)。M-K 突变 检验 表明 (图 
2b) ,整体 上 径流 呈现 由 多 变 少 的 演变 ,UF 和 UBi 交 
点 位 于 1992 年 ,说明 径流 模 数 的 突变 年 份 发 生 在 
1992 年 。 
2.1.2 给 沙 模 数 年 际 变化 特征 ”1957 一 2021 年 关 川 
河上 游 年 均 输 沙 模 数 为 2654.60 te km? +a, eK Al 
沙 模 数 为 11081.788 t-km- a! (1973 年) ,最 小 输 沙 
模 数 为 1.33 tkm?*…a'(2016 年 )。 区 间 内 输 沙 模 数 
星 波动 减少 趋势 (图 3a)。1957 一 1973 年 输 沙 模 数 
呈 增 加 趋势 (图 3a) ;1973 一 2021 年 输 沙 模 数 呈 减少 
趋势 ,其 中 1999 一 2021 年 输 沙 模 数 呈 显著 下 降 趋势 
(P<0.05)。 同 时 由 图 3b 可 知 ,UF 和 UB 的 交点 位 于 
2003 年 ,说 明 输 沙 模 数 在 2003 年 发 生 突变 。 
2.2 降水 量 特征 
2.2.1 降水 量 年 际 变 化 ”1957 一 2021 年 关 川 河上 游 
面 平均 年 降水 量 为 413.0 mm, 最 大 年 降水 量 为 
715.6 mm(1967 年) ,最 小 值 为 268.4 mm (1997 年) 
(图 4a) ,两 者 相差 302.6 mm, 最 大 年 降水 量 是 最 小 
年 降水 量 的 1.73 倍 。 降 水 量 总 体 呈 不 显著 下 降 趋 
$&(P>0.05) 0 UF, Fl UB, 26 F 2014 4F All 2020 4F , 2 
明 降水 量 突变 年 为 2014 年 和 2020 年 (图 4b)。 
2.2.2 降水 量 周 期 性 

(1) 周期 性 

由 图 5a 小波 实 部 等 值 线 分 布 可 知 , 降 水 在 不 同 
时 间 尺 度 上 存在 不 同 周期 丰 枯 循环 交 蔡 的 现象 。 
降水 量 的 周期 变化 存在 3 个 时 间 尺 度 2~4 a、4~9 a 
和 9~18 a。 其 中 ,2~4 a 在 小 时 间 尺 度 上 的 降水 振荡 
剧烈 ,周期 变化 规律 较为 清晰 ,信号 表现 稳定 ;4~9 a 
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图 2 径流 模 数 年 际 分 布 (a) 和 M-K 检 验 (b) 
Fig. 2 Annual distribution of runoff modulus (a) and M-K test (b) 
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图 3 输 沙 模 数 年 际 分 布 (a) 和 M-K 检 验 (b) 


Fig.3 Annual distribution of sediment transport modulus (a) and M-K test (b) 


1768 T T 
(a) 一 一 年 输 沙 模 数 
年 均 输 沙 模 数 
一 100 | e 年 输 沙 模 数 线性 趋势 
g 5000 AMIL, i 
; | Vy Ue 
S00 1960 1970 1980 1990 2000 2010 2020 
年 份 
1000 
(a) 一 一 年 降水 量 
年 均 降 水 量 
800 ate 年 降水 线性 趋势 
g 600 | 
a nn 人 
加 Tap VY ead ca 
200 


? 1960 1970 1980 1990 2000 2010 2020 
年 份 


4 1 1 L 1 1 L 1 
1960 1970 1980 1990 2000 2010 2020 
年 份 


图 4 降水 量 年 际 变 化 (a) 和 M-K 检 验 (b) 


Fig.4 Interannual variation of precipitation (a) and M-K test (b) 


时 间 尺 度 上 周期 均匀 分 布 且 具有 显著 全 域 化 特征 ， 
存在 11 次 完整 的 丰 枯 交替 变化 。2021 年 后 正 值 实 
线 没有 完全 闭合 ,说 明 2021 年 后 降水 将 呈 起 伏 减 少 
趋势 ;9~18 a 时 间 尺 度 上 的 降水 振荡 剧烈 , 且 周 期 性 
变化 具有 局 部 性 , 仅 在 1957 一 1980 年 周期 交替 变化 
规律 较 清 晰 ,此 后 规律 消失 ,周期 振荡 趋 于 平缓 。 
9~18 a 时 间 尺 度 存在 两 次 全 域 性 丰 枯 交 蔡 变 化 ,两 
个 降水 增 大 振荡 重心 分 别 为 1967 年 和 1977 年 ,两 
个 降水 减 小 振荡 中 心 分 别 为 1960 年 和 1973 年 。 综 
上 所 述 ,研究 区 降水 变化 周期 随 不 同时 间 尺 度 发 生 
相应 的 变化 规律 ,表现 为 大 时 间 尺 度 与 小 时 间 尺 度 
互相 包含 。 

(2) 小 波 方 差 检 验 

由 小 波 方差 图 可 知 , 每 个 方差 峰值 对 应 的 周期 
尺度 代表 该 时 间 尺 度 的 主 周期 。 时 间 序 列 存在 5 个 
波峰 ,包括 3 个 显著 和 两 个 不 显著 峰值 ,不 显著 峰值 
可 忽略 不 计 , 其 中 显著 峰值 依次 对 应 3 a .6 a、13 a, 


时 间 尺 度 6a 对 应 峰值 最 大 ,表明 6 a 左右 周期 振幅 
最 强 ,同时 说 明 4~9 a 时 间 尺 度 是 降水 量 序列 变化 
的 第 一 主 周期 ;时 间 尺 度 3 a 对 应 第 二 峰值 ,表明 2~ 
4 a 时 间 尺 度 周期 振荡 次 之 ,是 第 二 主 周期 ;13 a 对 
应 第 三 峰值 ,周期 振荡 最 弱 , 其 效应 可 忽略 不 计 。 
结合 小 波 实 部 等 值 线 图 可 以 看 出 ,在 6a 对 应 的 4~9 
a 时 间 尺 度 上 ,2021 年 小 波 系数 实 部 等 值 线 为 负 , 表 
明 2021 年 降水 均 处 于 下 降 期 。 由 于 6a 作 为 第 一 主 
周期 的 振荡 强度 大 于 第 二 及 其 他 峰值 ,因此 ,未 来 
短 时 期 内 降雨 呈现 减少 的 概率 较 大 (图 5b)。 
2.3 水 土 保持 措施 特征 
2.3.1 水 土 保 持 措施 量 和 径流 模 数 、 输 沙 模 数 变化 
随 着 水 土 保持 措施 和 工程 项 目的 不 断 实施 , 关 川 
河上 游 各 项 水 土 保持 措施 数量 不 断 增加 ,水 土 保持 
治理 程度 逐年 提高 。 按 照 研 究 区 水 土 保持 措施 实 
施 力 度 的 大 小 ,大 体 分 为 三 个 阶段 ,分 别 是 1957 一 
1980 年 .1980 一 2000 年 和 2000 一 2021 年 (图 6a)。 
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图 5 年 降水 量 小 波 系数 实 部 等 值 线 (a) 和 人 小波 方差 (b) 


Fig. 5 Real part contour line (a) and wavelet variance (b) of annual precipitation wavelet coefficients 


1957—1980 年 水 土 保持 措施 以 梯田 和 造林 为 6560 hm?。 随 着 水 土 保 持 措施 面积 的 逐渐 增加 , 径 


FE ,治理 方法 分 散 单一 ,1980 年 水 十 保持 措施 总 面 
积 达 到 21249.24 hm’, 治理 程度 达 12.96% ,年 治理 率 
0.56% ,主要 水 土 保持 措施 为 梯田 ,其 面积 达到 
14961.47 hm , 占 水 土 保 持 总 面积 的 70.419%。 

1981—2000 年 为 综合 治理 阶段 ,国家 重点 投入 
水 土 保持 治理 措施 ,实施 了 国家 级 、 省 级 小 流域 综 
合 治理 . 坡 耕 地 建设 . 洲 地 坝 建 设 工程 ,治理 程度 显 
著 提 高 ,2000 年 水 土 保持 面积 达 92477.833 hm’, iA 
理 程 度 达 56.39% ,年 治理 率 2.17% ,其 中 梯田 面积 达 
42060.01 hm ,造林 面积 达 32795.55 hm ,种 草 总 面 
积 达 16202.26 hm’. 

2001—2021 年 为 平稳 提高 阶段 ,在原 有 基础 上 
开展 了 大 规模 退耕 还 林 ( 草 ) 工 程 ,2021 年 水 土 保持 


ee ita 图 6b)。 


影响 水 沙 变化 的 主要 因素 
2.4.1 降水 量 与 径流 模 数 、 输 沙 模 数 的 关系 ”影响 
径流 模 数 输 沙 模 数 的 主要 因素 是 降水 和 以 水 土 保 


持 措 施 为 主 的 人 类 活动 。 按 照 时 间 周 期 将 降水 量 
与 径流 模 数 、. 输 沙 模 数 的 时 间 序 列 划分 为 3 个 时 间 
段 ,第 I 阶段 1957 一 1980 年 ,第 工 阶 段 1981 一 2000 
年 和 第 了 阶段 2001 一 2021 年 。 不 同时 有 段 降水 量 、 径 
流 模 数 和 输 沙 模 数 平均 值 见 表 1。 由 回归 分 析 可 知 
( 表 2) ,在 第 工 阶段 径流 模 数 、. 输 沙 模 数 对 降水 量变 

化 响应 强烈 ,趋势 协同 性 强 , 相 关系 数 较 高 ,降水 量 
和 径流 模 数 . 输 沙 模 数 相关 性 均 达 到 显著 水 平 (P< 
0.01) ,降水 量 对 径流 模 数 . 输 沙 模 数 的 决定 系数 尼 分 


别 为 0.471.0.302, 第 工 阶段 径流 模 数 和 输 沙 模 数 平均 
值 分 别 为 19551.10 ma.km- .ar 和 4558.80tkm-.ar 
第 卫 阶 段 径流 模 数 输 沙 模 数 对 降水 量变 化 响应 同 


面积 达 13.05 X10‘ hm’, TERRIS 79.56% ,年 治理 
率 1.10% ,其 中 梯田 面积 达 60963.01 hm’ ,造林 面积 
达 43130 hm ,种 草 总 面积 达 19835 hm? , 封 育 面积 达 
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图 6 各 项 水 土 保持 措施 数量 变化 (a) 和 径流 模 数 、. 输 沙 模 数 变 化 (b) 


Fig. 6 Changes in the quantity of various soil and water conservation measures (a) and changes in runoff modulus and sediment 


transport modulus (b) 
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表 1 不 同时 段 年 降水 量 \ 径 流 模 数 和 输 沙 模 数 平均 值 


Tab.1 Average annual precipitation, runoff modulus, and sediment transport modulus at different time periods 


特征 值 I (1957 一 1980 年 ) TI (1981—20004F ) 亚 (2001 一 2021 年 ) 
降水 量 /mm 424.1 418.5 389.9 
径流 模 数 /ma km”. a) 19551.10 10465.21 2153.30 
输 沙 模 数 /(t*km?.a) 4558.80 2896.92 247.64 


表 2 不 同时 段 年 降水 量 与 径流 模 数 、 输 沙 模 数 相 关 关 系 


Tab. 2 The correlation between annual precipitation and runoff modulus and sediment transport 


modulus at different time periods 


相关 关系 时 段 (年 份 ) 模型 R R F Sig. 

降水 量 (x)- 径 流 模 数 (y,) I (1957 一 1980 年 ) yı = 64.69x, —7885.79 0.703 0.471 21.479 0.000" 
TI (1981—20004F ) yı = 46.03x, - 8798.39 0.674 0.454 14.976 0.001" 

Il (2001—2021 4) yı=2.13x, + 1306.18 0.108 0.012 0.224 0.641 
降水 量 (x,)- 输 沙 模 数 (y,) I (1957 一 1980 年 ) y=14.59x, — 1625.89 0.549 0.302 9.511 0.005™ 
TI (1981—20004F ) Y2 = 12.59%, - 2374.10 0.596 0.356 9.933 0.006" 

TI (2001—2021 4F) Y2 = 1.16x, - 204.98 0.296 0.088 1.825 0.193 
径流 模 数 (y)- 输 沙 模 数 I (1957 一 1980 年 ) ya = 0.25y, 一 405.99 0.881 0.775 75.951 0.000" 
(y2) Il (1981—2000 4) y2=0.29y, — 148.01 0.941 0.885 138.919 0.000” 
TI (2001—2021 4 ) y,=0.14y,- 57.51 0.718 0.515 20.165 0.000" 


方差 检验 值 ;Sig. 为 显著 性 检验 值 。** 为 P<0.01。 下 同 。 


样 强 烈 ,趋势 协同 性 强 ,但 相 比 于 第 工 阶 段 相关 系 
数 较 低 ,降水 量 和 径流 模 数 、 输 沙 模 数 均 达 到 显著 
水 平 (P<0.01), 降 水量 对 径流 模 数 输 沙 模 数 的 决定 
系数 记分 别 为 0.454、0.356, 第 工 阶 段 径流 模 数 和 输 
沙 模 数 平均 值 分 别 为 10465.21 m? .km .a 1! 和 
2896.92 tknm2.a753; 第 亚 阶 段 径流 模 数 、. 输 沙 模 数 对 
降水 量变 化 响应 迅速 减弱 ,趋势 协同 异化 ,降水 量 
与 径流 模 数 . 输 沙 模 数 相关 性 均 未 达到 显著 水 平 (P> 
0.05) ,降水 量 对 径流 模 数 、 输 沙 模 数 的 决定 系数 尼 分 
别 为 0.012.0.088 ,第 亚 阶 段 径 流 模 数 和 输 沙 模 数 平 
均值 分 别 为 2153.30 ms .km-.a 和 247.64 tekm”. as 
结果 表明 ,在 1957—2021 年 关 川 河上 游 流域 降水 对 
径流 模 数 和 输 沙 模 数 的 影响 随时 间 延 长 逐渐 降低 ; 


注 :x 为 降水 量 (mm) ;x; 为 水 土 保持 措施 面积 (km?) ;yi 为 径流 模 数 (mi + ke a!) ;y; 为 输 沙 模 数 (t.km*…a);R 为 相关 系数 ; 尺 为 贡献 系数 ;为 


径流 模 数 与 输 沙 模 数 相关 关系 协同 性 强 ,径流 模 数 
对 输 沙 模 数 的 决定 系数 尼 分 别 为 0.775 , 0.885 
0.515 , 均 达 到 极 显 著 水 平 (P<0.01)。 

2.4.2 水 土 保 持 措 施 面积 与 径流 模 数 、 输 沙 模 数 的 
关系 ”在 第 [阶段 径流 模 数 、 输 沙 模 数 对 水 土 保持 
措施 面积 变化 响应 都 不 强烈 ,未 达到 显著 水 平 (P> 
0.05) ,水 土 保持 措施 面积 对 径流 模 数 、 输 沙 模 数 的 
决定 系数 尼 分 别 为 0.041.0.026( 表 3); 第 开 阶 段 径 
流 模 数 和 输 沙 模 数 相对 于 第 T 阶段 对 水 土 保 持 措 
施 面积 变化 响应 较 强 ,但 二 者 均 未 达到 显著 水 平 
(P>0.05) ,水 土 保持 措施 面积 对 径流 模 数 输 沙 模 数 
的 决定 系数 记分 别 为 0.158、0.033; 第 五 阶段 径流 模 
数 、 输 沙 模 数 对 水 土 保持 措施 面积 变化 响应 强烈 ， 


RI 不 同时 段 水 土 保 持 措施 量 和 径流 模 数 、 输 沙 模 数 相关 关系 


Tab. 3 The correlation between the amount of soil and water conservation measures and the runoff modulus and sediment 


transport modulus at different time periods 


相关 关系 时 段 (年 份 ) 模型 R R F Sig. 
水 土 保持 措施 面积 (x;)- 径 I (1957 一 1980 年 ) y, =-27.95x, +21337.21 0.202 0.041 0.939 0.343 
流 模 数 (7 ) [(1981 一 2000 年 ) y,=-10.29x, + 15914.73 0.398 0.158 3.385 0.082 
TI (2001—2021 4F ) yı =—7.40x, + 10519.78 0.673 0.454 15.771 0.001” 
水 土 保持 措施 面积 (x,)- 输 I (1957 一 1980 年 ) y= —6.38x, + 4966.52 0.160 0.026 0.579 0.455 
沙 模 数 (>) Il (1981—20002F:) y, =-1.445x, + 3661.78 0.181 0.033 0.067 0.446 
亚 (2001 一 2021 年 ) Y2 =-1.57x, 一 2024.07 0.717 0.514 20.059 0.000" 
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均 达 到 极 显 著 水 平 (P<0.01) ,水 土 保持 措施 面积 对 
径流 模 数 、. 输 沙 模 数 的 决定 系数 尼 分 别 为 0.454、 
0.514。 径 流 模 数 和 输 沙 模 数 平均 值 分 别 由 1957 一 
1980 年 的 19551.10 m°- km? +a” F1 4558.80 t-km?+a! 
减少 至 2001—2021 年 的 2153.30 mè- km~?» a Ail 
247.64t* km *…a!。 表 明 水 土 保持 措施 对 径流 模 数 
和 输 沙 模 数 的 影响 随时 间 延 长 逐渐 增 大 , 随 着 水 土 
保持 措施 面积 不 断 增加 ,水 土 流 失 得 到 有 效 控制 。 
2.5 年 降水 量 、 水 土 保持 措施 面积 对 年 径流 模 数 、 
年 输 沙 模 数 影响 路 径 

251 模型 初步 拟 合 检验 采用 结构 方程 模型 
(SEM ) 量 化 解释 变量 对 响应 变量 的 直接 影响 和 间接 
影响 ,对 年 径流 模 数 .年 输 沙 模 数 与 年 降水 量 7K 
保持 措施 面积 的 关系 进行 验证 分 析 , 初 步 拟 合 检验 
确定 数据 的 科学 性 .可靠 性 。 评 价 标准 为 :近似 卡 
方 值 (CMIN) 小 于 3, 拟 合 度 可 接受 ; 拟 合 优 度 指数 
(GFI) .相对 拟 合 指数 (NFT) .比较 拟 合 指数 (CEFT) 均 
大 于 0.8, 拟 合 度 很 好 ;近似 误差 均 方 根 (RMSEA) 小 
于 0.8, 拟 合 度 很 好 ( 表 4)。 各 项 指标 均 达 到 拟 合 标 
准 , 年 降水 量 、 年 径流 模 数 、 年 输 沙 模 数 水 土 保持 


的 间接 影响 显著 ,说 明 径 流 模 数 在 降水 量 水土 保 
持 措施 面积 和 输 沙 模 数 之 间 具 有 部 分 中 介 作 用 。 
模型 路 径 分 析 结 果 如 下 。 

关 川 河上 游 1957 一 2021 年 降水 量 与 径流 模 数 
的 关系 中 ,路径 系 数 的 标准 化 估计 值 为 0.42 ,表明 降 
水 量 每 增加 1 个 标准 单位 ,径流 模 数 增加 0.42 个 标准 
单位 ;降水 量 与 输 沙 模 数 的 关系 中 ,降水 量 对 输 沙 模 
数 的 直接 效应 不 显著 ;水 土 保持 措施 面积 与 径流 模 
数 的 关系 中 ,路 径 系数 的 标准 化 估计 值 为 -0.72, 表 明 
水 土 保 持 措 施 面积 每 增加 1 个 标准 单位 ,径流 模 数 
减少 0.72 个 标准 单位 ;水 土 保持 措施 面积 与 输 沙 模 
数 的 关系 中 ,水 土 保持 措施 面积 对 输 沙 模 数 的 直 
接 效应 不 显著 ;径流 模 数 与 输 沙 模 数 的 关系 中 ， 
路 径 系 数 的 标准 化 估计 值 为 0.99 ,表明 径流 模 数 
每 增加 1 个 标准 单位 , 输 沙 模 数 增加 0.99 个 标准 
单位 (图 7)。 

通过 模型 分 析 结 果 可 知 ( 图 7、 表 6) ,降水 量 与 
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$ 
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措施 面积 4 项 研究 数据 具有 较 好 的 拟 合 优 度 和 较 高 
的 可 靠 性 ,模型 适 配 良好 。 具 体 指标 结果 如 下 : 

2.5.2 水 沙 变化 路 径 ”路 径 分 析 结 果 表明 ( 表 5) , 径 
流 模 数 对 输 沙 模 数 起 到 决定 性 作用 ,同时 径流 模 数 
也 受到 降水 与 水 土 保 持 措施 的 影响 。 降 水 量 对 径 
流 模 数 的 直接 影响 显著 ,水 土 保持 措施 面积 对 输 沙 
模 数 的 直接 影响 显著 ,径流 模 数 对 输 沙 模 数 的 直接 
影响 显著 ,降水 量 与 水 土 保持 措施 面积 对 输 沙 模 数 


表 4 指标 初步 拟 合 检验 
Tab. 4 Preliminary fitting test of indicators 
拟 合 指标 ”CMIN GFI NFI CFI RMSEA 
仿 验 值 2.325 0.982 0.989 0.994 0.144 


表 5 路 径 关 系 检验 


Tab.S Path relationship verification 


变量 间 的 关系 路 径 系数 P 
降水 量 一 径流 模 数 0.417 0.000™ 
水 土 保持 措施 面积 一 径流 模 数 -0.720 0.000" 
径流 模 数 一 输 沙 模 数 0.993 0.000™ 
降水 量 一 输 沙 模 数 -0.039 0.491 
水 土 保持 措施 面积 一 输 沙 模 数 0.062 0.402 


注 :*** 表 示 P<0.001 时 有 极 显著 性 。 


, -0.72 | 
水 土 保持 措施 面积 0.06 


注 : 实 线 表 示 正 向 影响 ;虚线 表示 负 向 影响 ; 线 宽 表 示 影 响 大 小 。 
图 7 水 沙 变化 形成 的 结构 方程 模型 
Fig.7 The structural equation model for the formation of 


water and sediment variations 


水 土 保 持 措 施 均 对 径流 模 数 的 直接 影响 显著 ,对 输 
沙 模 数 的 直接 影响 不 显著 。 降 水 量 对 径流 模 数 和 
输 沙 模 数 的 影响 较 小 ,总 效应 分 别 为 0.42 和 0.38 ,其 
中 ,直接 效应 分 别 为 0.42 -0.04, 间 接 效应 分 别 为 0、 


表 6 效应 值 计 算 


Tab.6 Effect value calculation 


T > 
径流 模 数 输 沙 模 数 

降水 量 直接 效应 0.42 -0.04 
间接 效应 0 0.42 

总 效应 0.42 0.38 

水 土 保 持 措 施 面 积 ”直接 效应 -0.72 0.06 
间接 效应 0 -0.71 

总 效应 -0.72 -0.65 


注 : 直 接 影响 即 自 变 量 直 接 指向 因 变 量 的 路 径 系数 ,间接 影响 即 从 
自 变 量 指向 因 变 量 路 径 系 数 的 乘积 ,总 影响 即 直接 影响 和 间接 影 
响 的 总 和 。 


1772 干 旱 


40 卷 


0.42; 水 土 保持 措施 面积 对 径流 模 数 和 输 沙 模 数 的 
影响 最 大 ,总 效应 分 别 为 -0.72 和 -0.65, 其 中 ,直接 
效应 分 别 为 -0.72 .0.06 ,间接 效应 分 别 为 0.-0.71。 
主要 得 益 于 水 土 保持 综合 治理 ,退耕 还 林 还 草 等 生 
态 工程 的 实施 使 水 土 流失 状况 减轻 。 在 降水 量 减 
少 趋 势 不 显著 的 情况 下 ,径流 模 数 和 输 沙 模 数 呈现 
显著 减少 趋势 。 在 水 土 保持 措施 .降水 量 径流 模 
数 三 种 因素 中 径流 模 数 对 输 沙 模 数 起 到 决定 性 作 
用 ,由 此 可 推断 出 水 土 保持 措施 是 影响 输 沙 模 数 的 
最 重要 因素 ,通过 其 对 径流 模 数 的 直接 影响 间接 决 
定 输 沙 模 数 的 降低 。 


3 讨论 


关 川 河上 游 流域 降水 量 径流 模 数 和 输 沙 模 数 
变化 趋势 相同 ,但 降水 量 减 少 趋势 没有 达到 显著 水 
平 ;而 径流 模 数 、. 输 沙 模 数 均 呈 极 显 著 。 刘 强 汪 ”、 
张 富 2、 感 菲 写 等 对 不 同 流域 水 沙 变化 及 其 驱动 因 
素 进行 研究 ,结果 表明 水 沙 关 系 受 水 土 保 持 措 施 影 
响 变 化 显著 ,水 土 保持 措施 等 人 类 活动 可 以 有 效 降 
低 径 流 答 沙 量 。 王 进 呈 研究 发 现 祖 厉 河 流域 降水 
量 和 人 类 活动 对 径流 模 数 的 贡献 率 分 别 为 28.7% 和 
71.3% ,说 明 人 类 活动 是 重要 影响 因素 ,降水 量 是 次 
要 因素 ,这 与 本 文 的 研究 结果 一 致 。 本 研究 中 降水 
量 对 径流 模 数 的 决定 系数 尼 由 0.471 降 为 0.012 ,水 
土 保持 措施 面积 对 径流 模 数 的 决定 系数 让 由 0.041 
增 为 0.454; 降水 量 对 输 沙 模 数 的 决定 系数 尼 由 
0.302 降 为 0.088 ,水 土 保持 措施 面积 对 输 沙 模 数 的 
决定 系数 尼 由 0.026 增 为 0.514, 表 明 随 时 间 推 移 , 水 
沙 变化 的 主要 影响 因素 由 降水 量变 为 水 土 保持 措 
施 。 流 域内 水 土 保持 治理 程度 逐年 提高 ,由 1980 年 
的 12.96% 增 至 2021 年 的 79.56% ,有 效 降 低 了 径流 
模 数 和 输 沙 模 数 。 路 径 分 析 结 果 表 明 ,降水 和 水 土 
保持 措施 对 径流 模 数 . 输 沙 模 数 影响 的 路 径 是 通过 
直接 影响 径流 模 数 的 变化 ,进而 间接 影响 输 沙 模 数 
的 变化 。 从 总 效应 看 ,水土 保持 措施 面积 对 径流 模 
数 和 输 沙 模 数 总 效应 分 别 为 -0.72 .-0.65 ,降水 量 对 
径流 模 数 和 输 沙 模 数 的 总 效应 则 分 别 为 0.42 0.38, 
因而 水 土 保 持 措 施 是 影响 水 沙 变化 的 主要 因素 ,而 
降水 只 能 作为 次 要 因素 。 

气候 和 人 类 活动 引发 的 流域 水 沙 变化 一 直 是 
学 术 界 研究 的 热点 问题 。 在 众多 影响 水 沙 变化 的 


气候 影响 因子 和 人 类 活动 影响 因子 中 ,本 研究 仅 对 
降水 量 水 土 保持 措施 与 产 流产 沙 之 间 的 关系 进行 
了 探讨 ,其 他 因素 (温度 .土壤 流域 坡度 .季节 性 气 
候 变 化 、 降 十 强度 等 自然 因素 以 及 生产 建设 活动 等 
人 类 活动 ) 也 会 对 水 沙 变 化 产生 重要 影响 ,在 文中 
没有 给 予 充 分 考虑 ,有 待 此 后 进一步 深化 研究 。 
4 结论 

(1) 关 川 河上 游 流 域 径流 模 数 和 输 沙 模 数 在 研 
究 时 段 内 变化 同步 性 强 , 均 呈 显著 减少 趋势 ,其 突 
变 点 分 别 在 1992 年 和 2003 年 ;降水 量 在 1957 一 
2021 年 基本 稳定 , 呈 不 显著 减少 趋势 ;水 土 保持 措 
施 面积 总 体 随 时 间 推 移 不 断 增加 。 

(2) 径流 模 数 输 沙 模 数 对 降水 量变 化 的 响应 
由 强 转 弱 , 对 水 土 保持 措施 面积 的 响应 由 弱 变 强 。 
降水 量 对 径流 模 数 输 沙 模 数 的 相关 系数 及 分 别 由 
0.703 、0.549 减 小 为 0.108、0.296; 水 土 保持 措施 面积 
对 径流 模 数 输 沙 模 数 的 相关 系数 分 为 由 0.202、 
0.160 增 加 为 0.673、0.717, 说 明 水 土 保持 措施 面积 
对 流域 产 流产 沙 影响 逐渐 增 大 。 

(3) 降水 量 和 水 土 保持 措施 面积 对 径流 模 数 存 
在 直接 影响 ,径流 模 数 对 输 沙 模 数 存 在 直接 影响 ， 
降水 量 和 水 土 保持 措施 通过 直接 影响 径流 模 数 间 
接 影响 输 沙 模 数 。 降 水 量 和 水 土 保持 措施 对 水 沙 
变化 影响 的 关键 路 径 是 “水 土 保持 措施 面积 一 径流 
模 数 ”, 其 次 是 “水 土 保持 措施 面积 一 径流 一 输 沙 模 
数 ”。 效 应 值 分 析 结 果 为 :水 土 保持 措施 面积 对 径 
流 模 数 (-0.72) 和 输 沙 模 数 (-0.65) 的 总 效应 大 于 降 
水 量 对 径流 模 数 (0.42) 和 输 沙 模 数 (0.38) 的 总 效 
应 。 表 明 水 土 保持 措施 面积 大 幅度 增加 是 水 沙 关 
系 变化 的 主要 因素 ,而 降水 为 次 要 因素 。 
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Abstract: In order to further explore the influence of precipitation and soil and water conservation measures on 
runoff sediment, and to gain an in-depth understanding of the driving factors of water and sand changes in small 
and medium-sized rivers of the Yellow River, M-K test, Morlet wavelet, linear regression and structural equation 
modeling were used to study the changes in precipitation and the area of soil and water conservation measures in 
the upper Guanchuan River from 1957-2021, as well as their the mechanism and process of action were studied. 
The results showed that annual precipitation showed a non-significant decreasing trend (P >0.05), runoff modulus 
and sand transport modulus showed a significant decreasing trend (P<0.05), and the area of soil and water 
conservation measures showed an increasing trend. The response of runoff modulus and sand transport modulus 
to precipitation got weakened, and the response to soil and water conservation measures was gradually 
strengthened; the key path of precipitation and soil and water conservation measures affecting the change of 
water and sand was “area of soil and water conservation measures—runoff modulus”, followed by “area of soil 
and water conservation measures—runoff—sand transport modulus”, then “precipitation—runoff modulus”, and 
finally “precipitation—runoff modulus—sand transport modulus”. The total effect of precipitation on water and 
sand change is 0.42 and 0.38 respectively, while the total effect of soil and water conservation measures on water 
and sand change is - 0.72 and - 0.65 respectively. The main factor affecting water and sand change in the 
watershed is the area of soil and water conservation measures, and the secondary factor is the amount of 
precipitation. 

Keywords: changes in water and sediment; soil and water conservation measures; structural equation model; 
upper Guanchuan River 


